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A  METHOD  FOR  ESTIMATING  THE  MECHANICAL  PROPERTIES  OF  A  SOLID 
MATERIAL  SUBJECTED  TO  SIGNIFICANT  COMPRESSIONAL  FORCES  — 
PART  I:  NUMERICAL  THEORETICAL  SOLUTION 
FOR  A  SINGLE  THICK  PLATE 

1.  INTRODUCTION 

The  mechanical  properties  of  materials  create  displacement  and  stress  fields  that  often 
contribute  significantly  to  the  static  and  dynamic  response  of  the  structures  in  which  they  are 
found.  One  important  characteristic  shared  by  most  elastomeric  materials  (especially  slab¬ 
shaped  plates)  is  the  change  that  occurs  in  their  mechanical  properties  when  the  elastomer  is 
subjected  to  large  compressional  or  tensile  forces.  Under  these  forces,  the  rigidity  of  the 
material  typically  becomes  larger  and  the  damping  factors  become  smaller.  A  thorough 
understanding  of  such  behavior  is  necessary  so  that  the  static  and  dynamic  responses  of  a 
material  can  be  correctly  included  in  mathematical  models,  as  well  as  properly  understood  in 
the  actual  physical  structure  itself. 

Some  of  the  numerous  methods  developed  to  determine  the  properties  of  various  materials, 
include  those  based  on  the  use  of  resonance1 4  and  transfer  function  data/’5'2  Several  parameter 
estimation  techniques  have  also  been  investigated  for  plates.8'11  Although  the  above  approaches 
do  not  allow  for  testing  under  significant  compressional  forces,  efforts  have  been  made  to 
measure  material  properties  under  large  pressures.12'14  In  such  research,  the  material  is  placed  in 
a  pressurized  setting  and  insonified,  after  which  its  response  is  measured.  However,  these 
procedures,  which  are  typically  conducted  under  extreme  atmospheric  pressure  in  the  laboratory, 
can  have  an  adverse  effect  on  instrumentation,  as  well  as  on  safety.  In  other  studies,15’16  a  mass- 
loaded,  long,  thin  rod  has  been  examined  with  respect  to  the  bar  wavespeed  and  corresponding 
Young’s  modulus  (shear  motion  is  not  addressed). 

This  report  describes  an  inverse  method  that  has  been  developed  to  measure  complex, 
frequency-dependent  dilatational  and  shear  wavenumbers  of  a  single  slab-shaped  material 
subjected  to  large  compressional  forces.  Based  on  thick  plate  theory,  the  linear  equations  of 
motion  of  the  system  are  first  derived  for  a  test  specimen  that  is  attached  to  a  shaker  at  the 
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bottom  and  a  mass  at  the  top.  A  typical  test  configuration  is  shown  in  figures  1  and  2,  where  the 
shaker  projects  mechanical  energy  onto  a  plate-shaped  material  that  is  mass  loaded.  Two 
transfer  function  measurements  are  obtained  by  vibrating  the  mass-loaded  material  in  both  the 
vertical  and  horizontal  directions.  Once  this  process  is  accomplished,  the  transfer  functions  are 
combined  to  yield  closed-form  values  of  the  dilatational  and  shear  wavenumbers  at  any  given 
test  frequency.  After  these  parameters  are  estimated,  calculations  can  be  made  for  the  complex, 
frequency-dependent  dilatational  and  shear  wavespeeds;  Young’s  and  shear  moduli;  and 
Poisson’s  ratio. 

The  above  method  is  intended  for  use  with  materials  that  are  to  be  placed  in  an  environment 
where  they  will  be  subjected  to  large  compressional  forces,  such  as  those  that  would  typically 
arise  in  submarines,  where  the  panels  that  coat  the  exterior  of  the  ship  are  exposed  to  a  wide 
range  of  hydrostatic  pressures. 


Figure  1.  Vertical  Motion  Test 
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Figure  2.  Horizontal  Motion  Test 


2.  SYSTEM  MODEL 


The  motion  of  the  test  specimen  shown  in  figures  1  and  2  is  governed  by  the  equation 


/N2u  +  (A  + ju)VV»u  ~  P~^r  ,  (1) 

where  X  and  p  are  the  complex  Lame  constants  (N/m2 ),  p  is  the  density  ( kg/m3 ),  t  is 
time  (s),  •  denotes  a  vector  dot  product,  and  u  is  the  Cartesian  coordinate  displacement  vector 
of  the  material.  The  coordinate  system  of  the  test  configuration  is  shown  in  figure  3.  Note 
that  using  this  orientation  results  in  b  =  0  and  a  having  a  value  less  than  zero.  The  thickness 
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Figure  3.  Coordinate  System  of  Test  Specimen  Used  in  Model 


of  specimen  h  is  a  positive  value.  Equation  (1)  is  manipulated  by  writing  the  displacement 
vector  u  as 


\ux(x,y,z,tj\ 


u  = 


uy{x,y,z,t)  \  , 

uz(x,y,z,t) 


(2) 


where  x  is  the  location  along  the  plate  (m),  y  is  the  location  into  the  plate  (m),  and  z  is  the 
location  normal  to  the  plate  (m),  as  shown  in  the  figure.  The  symbol  V  in  equation  ( 1 )  is  the 
gradient  vector  differential  operator  written  in  three-dimensional  Cartesian  coordinates  as 


V  = 


d  .  d  . 

—  lr  H - / 

<?x  *  dy 


y 


(3) 


with  ix  denoting  the  unit  vector  in  the  x-direction,  iy  denoting  the  unit  vector  in  the  y-direction, 
and  /2  denoting  the  unit  vector  in  the  2-direction;  V2  is  the  three-dimensional  Laplace  operator 
operating  on  vector  u  as 


V2u  =  VVi  +  V2mv/'  +  V2u.i. 

x  a  }  y  4.  * 


(4) 
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and  operating  on  scalar  u  as 


V^u  =  V  •  Vw 

x,ytz  v  ™x,y,z 


dx 2  2  a - 2 


dy  dz 


The  term  V  *u  is  called  the  divergence  and  is  equal  to 


dux  duy  du. 


V  •  u  = — -  + — -  +  - 


dx  dy  dz 


(5) 


(6) 


The  displacement  vector  u  is  written  as 

u=V^+Vx^,  ^ 

wrhere  ^  is  a  dilatational  scalar  potential,  x  denotes  a  vector  crossproduct,  and  y/  is  an  equi- 
voluminal  vector  potential  expressed  as 


y/x{x,y,z,t) 
V=\y/y(x,y,z,t) 
xyz{x,y,z,t)\ 


(8) 


The  problem  is  formulated  as  a  two-dimensional  system;  thus  y  =  0,  uy{x,y,z,t)  s  0,  and 

d{  )!dy  =  0 .  Expanding  equation  (7)  and  breaking  the  displacement  vector  into  its  individual 
nonzero  terms  yields 

(9) 


(10) 

Equations  (9)  and  (10)  are  next  inserted  into  equation  (1),  which  results  in 


dx  dz 

and 

dz  dx 
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(11) 


and 


c] :V2y/  (x,z,t)- 


St2 


(12) 


where  equation  (1 1)  corresponds  to  the  dilatational  component  and  equation  (12)  corresponds  to 
the  shear  component  of  the  displacement  field.  Correspondingly,  the  constants  ctl  and  cs  are  the 
complex  dilatational  and  shear  wavespeeds,  respectively,  and  are  determined  by 


c 


d 


X  +  2/u 
P 


(13) 


and 


cc 


(14) 


The  relationship  of  the  Lame  constants  to  the  Young’s  and  shear  moduli  is  shown  as 


A  = 


Ev 

(l  +  o)(l-2v) 


(15) 


and 


P  =  G  = 


E 

2(l  +  o)  ’ 


(16) 


where  E  is  the  complex  Young’s  modulus  ( N/m2 ),  G  is  the  complex  shear  modulus  ( N/m2 ),  and 
v  is  the  Poisson's  ratio  of  the  material  (dimensionless). 

The  conditions  of  infinite  length  and  steady-state  response  are  now  imposed,  allowing  the 
scalar  and  vector  potential  to  be  written  as 

<j>{x,z,t)  =  ®(z)exp(i£x)exp(iroO  (17) 
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and 


¥y{x,z,t)  =  'F(z)exp(ifoc)exp(i©/)  ,  (18) 

where  i  is  the  square  root  of— 1,  ©is  frequency  (rad/s),  and  k  is  wavenumber  with  respect  to  the 
x-axis  (rad/m).  Inserting  equation  (17)  into  equation  (1 1)  yields 


The  solution  to  equation  (19)  is 

$0)  =  A(k,  ©)exp(i az)  +  B(k,  eo)  exp(-iaz)  , 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 
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and  the  solution  to  equation  (22)  is 


¥(z)  =  C(k,  <y)exp(i/?z )  +  D(k,  <y)exp(-i/?z)  ,  (26) 

where  A(k,co),  B{k,co),  C(k,co ),  and  D(k,co )  are  wave  response  coefficients  that  are  determined 
later  in  section  3.  The  displacements,  now  written  as  functions  of  the  unknown  constants  using 
the  expressions  in  equations  (9)  and  (10),  are 

u.(x,z,t )  =  U.(k,z,o))exp(\kx)exp(\o){) 

=  {\a[A(k,co)exp(\az)  -  B(k,co)Qxp(-\azj\  (27) 

+  ik[C(k,a>)exp(i/3z)  +  D(^«)exp(-i/?z)]}exp(i&x)exp(i&>0 


and 


ux(x,z,t )  =  Ux(k,z,  a)exp(ikx)exp(ia>t) 

=  \\k\A(k,co)exp(\az)+  B{k,co)exp(-\az'^  (28) 

-  i/?[C(A:,iy)exp(i/?z)  -  £>(&,&>)  exp(-i/?z)]  }exp(iAx)exp(ift;t)  . 

In  the  next  step,  specific  boundary  conditions  must  be  provided  to  obtain  the  individual  solutions 
for  vertical  motion  (section  3)  and  horizontal  motion  (section  4). 
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3.  TRANSFER  FUNCTION  FOR  VERTICAL  MOTION 


For  the  case  of  vertical  motion,  the  base  at  z  —  a  is  vibrated  vertically  with  a  shaker  (see 
figure  1).  Formulating  this  problem  requires  definition  of  the  four  boundary  conditions,  as 
shown  next. 

Because  a  rigid  mass  is  attached  to  the  material  and  the  particle  motion  is  vertical,  the 
tangential  (horizontal)  motion  at  the  top  of  the  plate  (z  =  b )  is  zero,  with  this  equation  written  as 

ux(x,b,t)  =  0  .  (29) 

The  normal  stress  at  the  top  of  the  specimen  is  equal  to  the  opposite  of  the  load  created  by  the 
mass  in  the  "-direction.  This  expression  is 

xzz  = a + 2fi)-uXx-bJ) + x — = -m—  hiiiM  m 

oz  dx  dt1' 

where  M  is  mass  per  unit  area  ( kg/m2 )  of  the  attached  mass.  The  tangential  motion  at  the 
bottom  of  the  plate  (z  =  a)  is  zero  and  is  shown  as 

ux(x,a,t)  =  0.  (3i) 

The  normal  motion  at  the  bottom  of  the  plate  (prescribed  as  a  system  input)  is  written  as 

uz  (x,  a,t)  =  U0  exp(icut)  .  (32) 

Assembling  equations  (1)— (32)  and  letting  b  =  0  yields  the  four-by-four  set  of  linear 
equations  that  model  the  system  as 


Ax  =  b 


where  the  entries  of  equation  (33)  are 


An  =  ik 


^12  —  5 


Au  =  -\fi 


m 

1 

II 

"3- 

(37) 

A2l  =  -a2A  - 2 a2/u  -  Ak2  -  i Mara  , 

(38) 

A22  = -a2A -2a2 /u- Ak2 +\Mco2a  , 

(39) 

A23  =  -2 k/5/j.  -  \Ma)2k  , 

(40) 

A24  =2kp/u-\Mco2k  , 

(41) 

Am  =  Auexp(iaa)  , 

(42) 

A32  =  Au  exp(-i  a  a)  , 

(43) 

A  33  =  A  l3exp(i/}ci)  , 

(44) 

A34  =  -A,  3  exp(-i fid)  , 

(45) 

A4I  =  iaexp(iaa)  , 

(46) 

A42  =  -iaexp(-iaa)  , 

(47) 

A43  =  ikexp(ij3a)  , 

(48) 

A44  =  \kexp(-\pa)  , 

(49) 

*11  =  A(k,a>)  , 

(50) 

*21  =  B(k,eo)  , 

(51) 

*31  =C(k,a>)  , 

(52) 

*4i  =  D(k,co)  , 

(53) 

S** 

1! 

o 

(54) 

Z>21  —  0  , 

(55) 

i3]  =  0  , 

(56) 

K  =  uQ . 

(57) 

From  equations  (34)-(57),  the  solution  to  the  constants  A{k,co),  B(k,oj),  C{k,co),  and  D{k,co) 
can  be  calculated  at  each  specific  wavenumber  and  frequency  from 

x  =  A  'b  .  (58) 
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Noting  that  k  -  0  for  vertical  motion  and  using  the  coefficients  from  equation  (58)  results  in  the 
transfer  function  between  the  vertical  base  displacement  and  the  vertical  mass  displacement 
being  written  as 


Tm= 


1  UMKa) 

Rx(o>)  U0 


cos(kdh) 


\P  J 


kd  sin (kdh) 


(59) 


where  Tx  (co)  or  Rx  ( co )  corresponds  to  the  frequency-domain  data  from  the  vertical  motion 

experiment,  which  are  typically  obtained  by  applying  a  Fourier  transform  to  raw  time-domain 
data  collected  with  accelerometers  or  laser  velocimeters. 


4.  TRANSFER  FUNCTION  FOR  HORIZONTAL  MOTION 

For  the  case  of  horizontal  motion,  the  base  at  z  =  a  is  vibrated  horizontally  with  a  shaker  (see 

fiourc  2).  Formulating  this  problem  requires  definition  of  the  four  boundary  conditions,  as 
shown  next. 

Because  a  rigid  mass  is  attached  to  the  material  and  the  particle  motion  is  horizontal,  the 
shear  (tangential)  stress  at  the  top  of  the  plate  is  equal  to  the  opposite  of  the  load  created  by  the 
mass  in  the  x-direction.  This  expression  is 


Tzx(x,b,t)  =  n 


dujjx,b,t)  du_(x,b,t) 


dz 


dx 


=  —M 


d2ux  ( x,b,t ) 
~d? 


(60) 


where  M  is  the  mass  per  unit  area  (kg/m2)  of  the  attached  mass.  The  normal  motion  at  the  top 
of  the  plate  (z  =  b )  is  zero,  with  this  equation  written  as 

uz(x,b,t)  =  0.  (61) 

The  tangential  motion  at  the  bottom  of  the  plate  (z  =  a),  prescribed  as  a  system  input,  is  shown  as 


ux(x,a,t)  =  V0  exp(ia)t)  , 


(62) 
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and  the  normal  motion  at  the  bottom  of  the  plate  is  zero,  which  is  expressed  as 


u.(x,a,t)  -  0  .  (63) 

Assembling  equations  ( 1 )— (28)  and  (60)-(63)  and  letting  b  =  0  yields  the  four-by-four  set  of 
linear  equations  that  model  the  system  as 

Ax  =  b  ,  (64) 

where  the  entries  of  equation  (64)  are 


Au  -  -2juka  -io)2Mk  , 

(65) 

An  =2/jka  -  \a)2 Mk  , 

(66) 

Al3  =  nP2  -  /jk2  +  \a> 2Mfi  , 

(67) 

Al4  =  /up2  -  /uk2  -\co2 Mp  , 

(68) 

A2 1  =i  a  , 

(69) 

A  -  -A 

n12  —  /i2]  , 

(70) 

A  23  =  i  k  , 

(71) 

A24  =  A23  , 

(72) 

A3I  =A23exp(iaa)  , 

(73) 

A32  =  A2 3  exp(-iaa)  , 

(74) 

A33  =-ipexp(ipa)  , 

(75) 

A 34  =  iPexp(-ipa)  , 

(76) 

A4]  =  A2]  exp(iaa)  , 

(77) 

A 42  =  - A2 ,  exp(-iaa)  , 

(78) 

A43  =  A23  exp(i pa)  , 

(79) 

A  44  =  A23  exp(-i Pa)  , 

(80) 

xu  =  A(k,a>)  , 

(81) 

x2l  =B(k,co )  , 

(82) 

*31  =C(k,a))  , 

(83) 
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*41  =  Otk.  ru)  , 

*n=0, 

^21  =  0  , 


*3,  =  f'o  ■ 
and 

^41  =  0  . 


(84) 

(85) 

(86) 

(87) 

(88) 


Equations  (65)-(88)  allow  calculation  of  the  solution  to  the  constants  A{k.co) ,  B(k.o:>). 

C  (k.(o) .and  D(k.oj)  at  each  specific  wavenumber  and  frequency  from 

x  =  A  'h  .  (89) 

Noting  that  k  —  0  for  horizontal  motion  and  using  the  coefficients  from  equation  (89)  results  in 
the  transfer  function  between  the  horizontal  base  displacement  and  the  horizontal  mass 
displacement  being  written  as 


T2(co)  = 


R2(®) 


Ux(0,b,a>) 


cos  (ksh)  - 


f  ht 
P  j 


K  sin (ksh) 


(90) 


where  T2(co)  or  R2(co)  corresponds  to  the  frequency-domain  data  from  the  horizontal  motion 
experiment. 


5.  INVERSE  SOLUTION  FOR  SYSTEM  WITHOUT  MASS 

Although  designed  for  the  testing  of  materials  subjected  to  compressional  forces,  the 
experiment  can  also  be  conducted  without  these  forces,  as  described  next.  This  simpler  version 
without  the  mass  has  an  inverse  solution  that  is  a  closed-form  expression  —  which  does  not 
occur  when  the  mass  is  added  to  the  experiment  and  corresponding  analysis.  When  mass  is  not 
present  in  the  experiment,  equations  (59)  and  (90)  become 

T3(a?)  =  — —  =  _ I _ 

R3(co)  U0  cos  (kdh) 
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and 


R^co)  V0  cos  (ksh) 


(92) 


respectively.  Inverting  equation  (91)  yields  the  value  of  kd  as  a  function  of  /?3 .  The  solution 
for  the  real  part  of  kd  is 


Re(^)  = 


— Arccos(s)  +  — 
2  h  2  h 


— Arc  cos(-s)  +  — 
[2h  2  h 


n  even 


n  odd 


(93) 


where 

5  =  [Re(i?3 )]2  +  [Im(/?3 )]2  -  V{[R e(R3)]2  +  [Im(R3)]2F  -  {2[Re(R3)]2  "  2[Im(/?3)]2  -  l}  , 

(94) 

n  is  a  non-negative  integer,  and  capital  A  denotes  the  principal  value  of  the  inverse  cosine 
function.  The  value  of  n  is  determined  from  the  function  s,  which  is  a  cosine  function  with 
respect  to  frequency.  At  zero  frequency,  n  is  0.  Every  time  5  cycles  through  rr  radians  (1 80°), 
n  is  increased  by  1 .  After  the  solution  to  the  real  part  of  kd  is  found,  the  solution  to  the 
imaginary  part  is  written  as 


1  Re(J?3) 

Im(^)  =  -log  1 


lm  (*,) 


h  0  cos[R c{kd)h]  sin[R e(kd)h] 


(95) 


The  inverse  solution  to  equation  (92)  is  the  same  as  that  for  equation  (91)  and  is  written  as 


Re(*s)  = 


1  A  /A  m7t 

—  Arccos(r)  + -  m  even 

2  h  2h 


1  .  .  .  mn  ,, 

— Arccos(-r)  + -  m  odd 

2  h  2  h 


(96) 
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where 


r  =  [Re(i?4 )]2  +  [Im(i?4  )]2  -  ^/{[Re(^4  )]2  +  [Im(R4  )]2  f  -  {2[Re(i?4}]2  -2[Im(R4)]2  -l}  , 

(97) 

m  is  a  non-negative  integer,  and  capital  A  denotes  the  principal  value  of  the  inverse  cosine 
function.  The  value  of  m  is  determined  from  the  function  r,  which  is  a  cosine  function  with 
respect  to  frequency.  At  zero  frequency,  m  is  0.  Every  time  r  cycles  through  n  radians  (1 80°), 
m  is  increased  by  1 .  After  the  solution  to  the  real  part  of  ks  is  found,  the  solution  to  the 
imaginary  part  is  written  as 


Im(^)  =  — Ioge 
h 


Re(i?4) 

cos[Re(£s)/i] 


Im  (*4)  1 

sin[Re(&s  )h]  J  ' 


(98) 


6.  INVERSE  SOLUTION  FOR  SYSTEM  WITH  MASS 

Next  solved  is  the  inverse  problem  for  vertical  and  horizontal  motion  with  a  mass  attached  to 
the  plate.  This  procedure  involves  use  of  the  experimental  data  and  equations  (59)  and  (90)  to 
estimate  the  dilatational  and  shear  wavenumbers,  respectively.  Equations  (59)  and  (90)  can  be 
rewritten  as 

(99) 


(100) 

respectively,  where  the  problem  now  becomes  finding  the  zeros  of  the  right-hand  side  of 
equations  (99)  and  (100)  or,  in  the  presence  of  actual  data  that  contain  noise,  finding  the  relative 


f(kd)  =  0  =  cos(kdh)- 


\P) 


kd  s\n{kdh)  -  R{ 


and 


f(ks )  =  0  =  cos(ksh)  — 


f— 1 


ks  sin (ksh)  -  R 


2  5 
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minima  of  the  right-hand  side  of  equations  (99)  and  (100)  and  then  determining  which  relative 
minimum  corresponds  to  dilatational  and  shear  wave  propagation  and  which  relative  minima  are 
extraneous.  Because  equations  (99)  and  (100)  have  a  number  of  relative  minima,  zero-finding 
algorithms  are  not  applied  to  this  function,  as  they  typically  do  not  find  all  the  minima  locations 
and  are  highly  dependent  on  initial  starting  locations.  The  best  method  for  finding  all  minima 
locations  is  to  plot  the  absolute  value  of  the  right-hand  side  of  equations  (99)  and  (100)  as 
surfaces,  with  the  real  part  of  the  wavenumber  ( kd  or  ks )  on  one  axis  and  the  imaginary  part  of 
the  wavenumber  ( kd  or  ks )  on  the  other.  At  the  lower  frequencies,  the  minimum  farthest  to  the 
left  will  correspond  to  dilatational  or  shear  wave  propagation.  As  the  frequency  increases, 
extraneous  minima  will  appear  to  the  left  of  the  minimum  that  corresponds  to  dilatational  or 
shear  wave  propagation;  however,  the  wave  propagation  minimum  will  always  be  close  to  the 
previous  test  frequency  wave  propagation  minimum  provided  that  the  frequency  increments  are 
relatively  small. 

Also,  because  the  real  part  of  the  both  the  dilatational  and  shear  wavenumbers  are  mono- 
tonically  increasing  functions  with  respect  to  frequency,  each  increase  in  frequency  will  require 
that  the  real  part  of  the  new  wavenumber  that  is  estimated  be  greater  than  the  real  part  of  the  old 
wavenumber  that  was  previously  estimated.  Sometimes  referred  to  as  a  grid  method,  this  process 
is  illustrated  in  section  8. 
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7.  DETERMINATION  OF  PROPERTIES  FROM  WAVENUMBERS 


The  material  properties  can  be  calculated  from  the  wavenumbers.  First,  the  dilatational  and 
shear  wavespeeds  are  determined  from 


c 


d  ~ 


(101) 


and 


(102) 


respectively.  Next,  the  Lame  constants  are  calculated  from  equations  (13)  and  (14),  now  written 
as 


P=pc] 

and 

X=pc]  -2 pc]  . 

Poisson’s  ratio  is  then  determined  from 

o  = - , 

2(p  +  X) 

Young’s  modulus  from 

E=  2p^£  +  M) 

2{p  +  X)  ’ 

and  the  shear  modulus  directly  from 
G  =  p  . 


(103) 


(104) 


(105) 


(106) 


(107) 
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8.  NUMERICAL  SIMULATION 


The  above  measurement  method  can  be  simulated  by  means  of  a  numerical  example,  with  the 
values  for  the  soft  rubberlike  material  properties  of  the  test  specimen  expressed  as  follows: 
uncompressed  Young  s  modulus  E  of  [(le8  - i2e7)  +  ( 5e3  /~  —  i3e2 / ) |  in  N/m2  (where  /  is 
frequency  in  Hertz),  Poisson’s  ratio  v  of  0.40  (dimensionless),  density  p  of  1200  kg/m3 ,  and 
thickness  h  of  0.1  m.  The  top  mass  is  a  0,0254-m  (1-inch)  steel  plate  that  has  a  mass  per  unit 
area  value  M  of  199  kg/m".  With  this  mass  value,  the  compressed  Young’s  modulus  of  the  test 
specimen  is  1.5  times  the  uncompressed  value.  The  system  is  first  analyzed  without  the  mass, 
using  the  closed-form  solution  method  developed  in  section  5.  Next,  the  mass  is  added  to  the 
plate  and  the  grid  method  developed  in  section  6  is  applied  to  the  simulated  data. 

The  simulations  are  conducted  with  three  different  amounts  of  additive  noise:  no  additive 
noise,  2%  additive  noise,  and  4%  additive  noise.  Additive  noise  is  included  in  the  transfer 
function  with  the  equation 

Te  {(o)  =  T(cd)  +  e{Re[T (a)]cra  +  iIm[T (co)]ab }  ,  (108) 

where  e  is  the  amount  of  additive  noise  added  to  the  transfer  function  and  cra  and  ab  are 

random  numbers  with  zero  mean  and  a  variance  of  one.  The  value  e  is  also  called  the  transfer 
function  error  value,  and  it  represents  the  deviation  of  the  transfer  function  from  perfect  data. 

Figure  4  plots  the  transfer  function  of  the  system  for  vertical  motion  without  the  mass  and 
corresponds  to  equation  (59)  with  M=  0.  Figure  5  illustrates  the  transfer  function  of  the  system 
for  horizontal  motion  without  the  mass  and  corresponds  to  equation  (90)  with  M=  0.  In  both 
figures,  the  top  plot  is  the  magnitude  in  decibels  and  the  bottom  plot  is  the  phase  angle  in 
degrees. 

Figure  6  is  a  plot  of  the  function  s  versus  frequency  and  corresponds  to  equation  (94).  Figure 
7  is  a  plot  of  the  function  r  versus  frequency  and  corresponds  to  equation  (97). 

All  the  plots  in  figures  4  through  7  are  displayed  without  additive  noise. 
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Figure  4.  Transfer  Function  of  Vertical  Motion  Versus  Frequency 
for  System  Without  Mass 
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Figure  8  shows  actual  and  estimated  dilatational  wavespeed  versus  frequency  for  a  system 
without  mass  for  no  additive  noise,  2%  additive  noise,  and  4%  additive  noise.  This  plot 
corresponds  to  the  solution  determined  from  equations  (93)— (95).  Figure  9  shows  actual  and 
estimated  shear  wavespeed  versus  frequency  without  mass  for  no  additive  noise,  2%  additive 
noise,  and  4%  additive  noise.  This  plot  corresponds  to  the  solution  determined  from  equations 
(96)— (98). 

Figure  10  plots  actual  and  estimated  shear  modulus  versus  frequency  for  a  system  without 
mass  for  no  additive  noise,  2%  additive  noise,  and  4%  additive  noise  and  is  constructed  using 
equation  (103).  Figure  1 1  shows  actual  and  estimated  Young’s  modulus  versus  frequency  for  a 
system  without  mass  for  no  additive  noise,  2%  additive  noise,  and  4%  additive  noise  and  is 
constructed  using  equations  (103),  (104),  and  (106). 

Figure  12  illustrates  actual  and  estimated  real  Poisson’s  ratio  versus  frequency  for  a  system 
without  mass  for  no  additive  noise,  2%  additive  noise,  and  4%  additive  noise  and  was  plotted 
using  equation  (105).  Because  the  numerical  example  is  formulated  using  a  Poisson’s  ratio  that 
is  strictly  real,  no  imaginary  component  is  shown  in  this  plot.  Imaginary  values  of  Poisson’s 
ratio  are  possible,  however,  and  have  been  shown  to  theoretically  exist.17 

Figure  1 3  is  a  plot  of  the  transfer  function  of  the  system  with  mass  for  vertical  motion  versus 
frequency  and  corresponds  to  equation  (59)  with  M=  199  kg/m2 .  Figure  14  is  a  plot  of  the 
transfer  function  of  the  system  with  mass  for  horizontal  motion  versus  frequency  and  corre¬ 
sponds  to  equation  (90)  with  M=  199  kg/m" .  The  plots  in  both  figures  are  displayed  without 
additive  noise. 

Figure  15  shows  the  magnitude  of  the  surface  given  by  equation  (99)  versus  real  and 
imaginary  dilatational  wavenumbers  at  2000  Hz  with  an  additive  noise  value  of  0.04  for  a  system 
with  mass.  The  actual  value  of  the  dilatational  wavenumber  is  denoted  with  a  square  marker,  and 
the  estimated  value  (found  using  the  identification  of  a  local  minimum)  is  shown  on  the  plot  with 
a  circle  marker.  Figure  16  is  a  plot  of  the  magnitude  of  the  surface  given  by  equation  (100) 
versus  real  and  imaginary  shear  wavenumbers  at  2000  Hz  with  an  additive  noise  value  of  0.04  for 
a  system  with  mass.  The  actual  value  of  the  shear  wavenumber  is  denoted  with  a  square  marker, 
and  the  estimated  value  (found  using  the  identification  of  a  local  minimum)  is  shown  on  the  plot 


23 


Imaginary  Dilatational  Wavespeed  (m/s)  Real  QMatational  Wavespeed  (m/s) 


20 

0 

-20 

-40 

-60 

-80 

■100 

-120 

0  1000  2000  3000  4000  500 

Frequency (Hz) 

Figure  8.  Actual  and  Estimated  Dilatational  Wavespeed  Versus  Frequency 

for  System  Without  Mass 
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Figure  9.  Actual  and  Estimated  Shear  Wavespeed  Versus  Frequency 
for  System  Without  Mass 
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Figure  11.  Actual  and  Estimated  Young’s  Modulus  Versus  Frequency 

for  System  Without  Mass 
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Figure  12.  Actual  and  Estimated  Rea /  Poisson ’s  Ratio 
Versus  Frequency  for  System  Without  Mass 


Figure  14.  Transfer  Function  of  Horizontal  Motion  Versus  Frequency 

for  System  with  Mass 
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Figure  15.  Magnitude  of  Surface  Versus  Real  and  Imaginary  Dilatational 
Wavenumbers  at  2000  Hz  Showing  Estimated  and  Actual  Values 
for  System  with  Mass  and  0.04  Additive  Noise 
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Figure  16.  Magnitude  of  Surface  Versus  Real  and  Imaginary  Shear 
Wavenumbers  at  2000  Hz  Showing  Estimated  and  Actual  Values 
for  System  with  Mass  and  0.04  Additive  Noise 
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with  a  circle  marker.  These  figures  were  constructed  using  200  discrete  points  in  both  axis 
directions.  The  identification  process  was  repeated  every  1 00  Hz  with  a  frequency  spectrum  of 
50  to  5000  Hz. 

Figure  17  is  a  plot  of  actual  and  estimated  dilatational  wavespeed  versus  frequency  for  a 
system  with  mass  for  no  additive  noise,  2%  additive  noise,  and  4%  additive  noise.  Figure  18 
shows  the  actual  and  estimated  shear  wavespeed  versus  frequency  for  a  system  with  mass  for  no 
additive  noise,  2%  additive  noise,  and  4%  additive  noise.  Figure  19  plots  actual  and  estimated 
shear  modulus  versus  frequency  for  a  system  with  mass  for  no  additive  noise,  2%  additive  noise, 
and  4%  additive  noise.  The  step-shaped  results  in  the  real  part  of  figures  18  and  19  (upper  plot) 
are  due  to  the  discretization  size  of  the  surfaces  in  equations  (99)  and  (100).  A  discretization  size 
finer  than  200-by-200  points  would  produce  less  steplike  (smoother)  results.  Figure  20  is  a  plot 
of  actual  and  estimated  Young’s  modulus  versus  frequency  for  a  system  with  mass  for  no 
additive  noise,  2%  additive  noise,  and  4%  additive  noise.  Figure  21  illustrates  actual  and 
estimated  real  Poisson’s  ratio  versus  frequency  for  a  system  with  mass  for  no  additive  noise,  2% 
additive  noise,  and  4%  additive  noise. 

Table  1  presents  additive  noise  versus  parameter  estimation  error  for  the  system  without 
mass  and  table  2  for  the  system  with  mass.  The  parameter  estimation  error  was  determined  from 
the  equation 

Q  _  J_  V"1  Kc.vf  (^/I  )  ~  (  J  09) 

k«K)| 

where  0is  the  parameter  estimation  error,  kcxI  (con )  is  the  estimated  value  of  the  parameter  at  the 
nth  frequency  value,  rcacl(o)n)  is  the  actual  value  of  the  parameter  at  the  nth  frequency  value,  and 
N  is  the  total  number  of  frequencies  at  which  an  estimate  is  computed.  Because  the  routine  did 
not  produce  a  realistic  estimate  of  the  parameters  at  very  low  frequencies,  these  parameters  are 
not  included  in  the  tables.  All  the  parameters  are  estimated  with  very  small  error.  Based  on  the 
results  shown  in  each  table,  it  is  found  that  the  grid  method  produces  a  more  accurate  estimate  of 
the  parameters  than  does  the  numerical  solution  method  for  data  with  noise. 


32 


(S/UJ)  paadsaABM  iBUOijejeiia  |Bay  (s/iu)  paadsaABM  leuojjeiena  Ajbu!6biu| 


Imaginary  Shear  Wavespeed  (m/s)  Real  Shear  Wavespeed  (m/s) 


Frequency (Hz) 


Frequency (Hz) 


Figure  18.  Actual  and  Estimated  Shear  Wavespeed  Versus  Frequency 

for  System  with  Mass 
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Figure  21.  Actual  and  Estimated  Real  Poisson ’s  Ratio 
Versus  Frequency  for  System  with  Mass 


Table  1.  Additive  Noise  Versus  Parameter  Estimation  Error 
for  System  Without  Mass 


Parameter 

Parameter 

Additive  Noise 

Estimation  Error 

Table  2.  Additive  Noise  Versus  Parameter  Estimation  Error 
for  System  with  Mass 


Parameter 

Additive  Noise 

Parameter 

Estimation  Error 

cd 

0 

0.0023 

cd 

0.020 

0.0076 

cd 

0.040 

0.013 

cs 

0 

0.0024 

cs 

0.020 

0.0046 

cs 

0.040 

0.0071 

G 

0 

0.0047 

G 

0.020 

0.0092 

G 

0.040 

0.014 

E 

0 

0.0045 

E 

0.020 

0.0087 

E 

0.040 

0.013 

V 

0 

0.0012 

V 

0.020 

0.0030 

u 

0.040 

0.0051 
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9.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  documents  the  derivation  of  a  theoretical  method  for  estimating  the  mechanical 
properties  of  slab-shaped  materials  subjected  to  compressional  forces.  After  a  single  thick  plate 
is  vibrated  both  vertically  and  horizontally  and  the  transfer  functions  between  the  faces  of  the 
plate  are  measured,  estimations  of  dilatational  wavespeed,  shear  wavespeed.  Young’s  modulus, 
shear  modulus,  and  Poisson’s  ratio  can  be  accurately  obtained.  Further  investigation  showed  that 
this  technique  is  relatively  immune  to  noise  mechanisms  that  are  sometimes  present  in  the 
measurements. 

It  is  recommended  that  actual  measured  data  be  used  in  a  laboratory  setting  to  evaluate  the 
effectiveness  of  the  inverse  method. 


10.  REFERENCES 


1 .  D.  M.  Norris  Jr.,  and  W.  C.  Young,  “Complex  Modulus  Measurements  by  Longitudinal 
Vibration  Testing,”  Experimental  Mechanics,  vol.  10,  1970,  pp.  93-96, 

2.  W.  M.  Madigosky  and  G.  F.  Lee,  “Improved  Resonance  Technique  for  Materials 
Characterization,’  Journal  of  the  Acoustical  Society  of  America,  vol.  73,  no.  4, 1983. 
pp.  1374-1377. 

3.  S.  O.  Oyadiji  and  G.  R.  Tomlinson,  “Determination  of  the  Complex  Moduli  of  Viscoelastic 
Structural  Elements  by  Resonance  and  Non-Resonance  Methods,”  Journal  of  Sound  and 
Vibration,  vol.  101,  no.  3,  1985,  pp.  277-298. 

4.  S.  L.  Garrett,  “Resonant  Acoustic  Determination  of  Elastic  Moduli,”  Journal  of  the 
Acoustical  Society  of  America,  vol.  88,  no.  1,  1990,  pp.  210-220. 

5.  J.  L.  Buchanan,  “Numerical  Solution  for  the  Dynamic  Moduli  of  a  Viscoelastic  Bar,” 
Journal  of  the  Acoustical  Society  of  America,  vol.  81,  no.  6,  1987,  pp.  1775-1786. 


41 


6.  T.-K.  Ahn  and  K.-J.  Kim,  “Sensitivity  Analysis  for  Estimation  of  Complex  Modulus  of 
Viscoelastic  Materials  by  Non-Resonance  Method,”  Journal  of  Sound  and  Vibration , 
vol.  176,  no.  4,  1994,  pp.  543-561. 

7.  A.  J.  Hull,  “An  Inverse  Method  to  Measure  the  Axial  Modulus  of  Composite  Materials 
Under  Tension  ''Journal  of  Sound  and  Vibration ,  vol.  195,  no.  4,  1996,  pp.  545-551. 

8.  S.  I.  Roklin  and  W.  Wang,  “Measurements  of  Elastic  Constants  of  Very  Thin  Anisotropic 
Plates,”  Journal  of  the  Acoustical  Society  of  America,  vol.  94,  no.  5,  1993,  pp.  2721-2730. 

9.  V.  Marchand,  J.  Authesserre,  J.  Pouyet,  and  C.  Bacon,  “Determination  of  the  Elastic 
Constants  of  Materials,  in  the  Form  of  Plates,  by  a  Free  Vibration  Method,”  Journal  of 
Sound  and  Vibration,  vol.  194,  no.  4,  1996,  pp.  497-512. 

10.  F.  Moussu  and  M.  Nivoit,  “Determination  of  Elastic  Constants  of  Orthotropic  Plates  by  a 
Modal  Analysis/Method  of  Superposition,”  Journal  of  Sound  and  Vibration,  vol.  165, 
no.  1,  1993,  pp.  149-163. 

11.  A.  J.  Hull,  “A  Method  for  Estimating  the  Mechanical  Properties  of  a  Solid  Material 
Subjected  to  Insonification  -  Part  I:  Theory,”  NUWC-NPT  Technical  Report  1 1,404, 
Naval  Undersea  Warfare  Center  Division,  Newport,  RI,  6  December  2002. 

12.  J.  E.  Cole  III  and  K.  Martini,  “Low-Frequency  Impedance  Measuring  Device  for 
Underwater  Acoustic  Material,”  Journal  of  the  Acoustical  Society  of  America,  vol.  95, 
no.  5,  part  2,  1 994,  p.  293 1 . 

13.  S.  Beretti,  C.  Audoly,  and  P.  Guillaussier,  “Pressure  Effects  on  the  Dynamic  and  Static 
Properties  of  Perforated  Elastomers  Used  in  Underwater  Acoustics,”  Journal  of  the 
Acoustical  Society  of  America,  vol.  97,  no.  5,  part  2,  1995,  p.  3293. 

14.  V.  E.  Glazanov  and  A.  V.  Mikhailov,  “Calculation  of  the  Velocity  of  Sound  in  Rubber 
Containing  Cylindrical  Voids  with  Allowance  for  Dynamic  Correction,”  Acoustical 
Physics,  vol.  47,  no.  6,  2001,  pp.  671-676. 

15.  J.  D.  Dickens,  “Phase  Velocity  of  Rubber  Element  in  Vibration  Isolator  Under  Static 
Load,”  Journal  of  Sound  and  Vibration,  vol.  234,  no.  1,  2000,  pp.  21-42. 

42 


1 6.  J.  D.  Dickens,  “Dynamic  Model  of  Vibration  Isolator  Under  Static  Load  Journal  of  Sound 
and  Vibration,  vol.  236,  no,  2,  2000,  pp.  323-337. 

17.  T.  Pritz,  “Frequency  Dependencies  of  Complex  Moduli  and  Complex  Poisson’ s  Ratio  or 
Real  Solid  Materials,’  Journal  of  Sound  and  Vibration,  vol.  214,  no.  1,  1998,  pp.  83-104, 


43  (44  blank) 


INITIAL  DISTRIBUTION  LIST 


Addressee 

No.  of  Copies 

Office  of  Naval  Intelligence  (ONI  241  -  J.  Zilius,  T.  Morgan- 
ONI  263 -S.  Brown  (2))  g  ’ 

Office  of  Naval  Research  (ONR  32 1  -  R.  Elswick;  ONR  333  -  K.  Ng)  2 

Defense  Technical  Information  Center  2 

Center  for  Naval  Analyses 


